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Abstract — A damage detection method is presented for the 
identification and quantification of damage. The proposed 
method uses finite element method to extract modal 
parameters of cracked and intact cantilever beam. The damage 
is simulated by fracture mechanics concept by introducing 
cracked elements at different locations with predetermined 
magnitude of depth. The curvature response function, function 
of crack location and size, are approximated by means of 
polynomial surface fitting. The numerical data obtained is 
meshed using B-spline. The algorithm based on curvature, 
Wavelet Transform and surface fitting technique is proposed 
for damage detection. Cubic fit and quadratic fit are used for 
function generation and three dimensional plots. Number of 
numerical examples is presented to demonstrate the accuracy 
of the proposed methodology. 

Index Terms — Crack beam modeling, Response Curvature 
damage surface, Wavelet Transform, Response Curvature 
error damage surface. 

I. Introduction 

The Vibration Based Inspection Methods (VBI) are 
popular because of convenience of measurement and 
collection of modal parameters. The global nature of the 
parameters (such as natural frequencies) implies that tests 
can be conducted at virtually arbitrary points. A prior 
knowledge of the damage location is not necessary. VBI 
methods also detect Hidden damage which is normally 
difficult to access. This capability draws attention from 
mathematicians and engineers during last decade. The crack 
beam modeling is important for numerical simulation. 
Modeling cracked beam element using the linear facture 
mechanics concepts [ 1 ] [2] will help in modeling hairline crack. 
An analytical as well as experimental method is used for crack 
detection in cantilever beams [3] . To locate the crack, contours 
of the normalized frequency in terms of the normalized crack 
depth and location are plotted. The intersection of contours 
with the constant modal natural frequency planes is used to 
relate the crack location and depth. The transverse crack is 
also modeled in a beam as rotational spring having stiffness 
of negligible mass [4] . For the computational convenience 
lumped mass matrix instead of consistent mass matrix is 
considered. The applications of Wavelet Transform to detect 
crack-like damage in structures are demonstrated. [5] [6] [7] 
This paper presents a method for damage detection based 
on the Response Curvature Damage Surface. The surface is 
generated from the curvature of mode shapes for first five 
modes. Additional data points are obtained by fitting B-spline 
curve. Curvature damage response surface function, 
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curvature error damage response surface function and 
curvature resultant damage response functions are obtained 
from polynomial fitting of numerical data. The wavelet 
transform is used to analyze curvature mode shapes which 
locate the crack. The algorithm addresses the determination 
of crack depth using the damage response surface fitting 
polynomial. 

II. Crack beam element modeling 

According to the principle of Saint- Venant's, the stress 
field is affected only in the region adjacent to the crack. It is 
difficult to find an appropriate shape function to approximate 
the kinetic energy and elastic potential energy because of 
the discontinuity of deformation in crack element. The 
calculation of the additional stress energy due to crack is 
computed using the concept in fracture mechanics .The 
flexibility coefficient expressed by a stress intensity factor 
is derived by applying the Castigliano's theorem in the 
linear-elastic range. A crack model beam can be divided into 
elements and the behavior of the elements located to the 
right of the cracked element may be regarded as external forces 
applied to the cracked element, while the behavior of 
elements situated to its left as constraints . Thus, the 
flexibility matrix of a cracked element with constraints is 
calculated. Neglecting shear action, the strain energy of an 
undamaged element is 
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Wbere E is the elastic modulus, P and M are the shear and 
bending internal forces at the right node, / the moment of 
inertia of the undamaged element and L the length of the 
finite element. For a rectangular beam having width b and 
thickness h additional strain energy due to the crack can be 
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Where K ,K & K are stress intensity factors for opening 
type, sliding type and tearing type cracks, respectively. Ep=E 
for plane stress, Ep=E/(l- v 2 ) for plane strain, and a is the 
crack depth. Taking into account only bending, equation (2) 
becomes 
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Where K K ,K are stress intensity factors for opening- 
type and sliding-type cracks due to M and P respectively 

K M = (6M I bk 1 ) 4^F } (s) 
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F(s) and F (s) are function of the ratio s between the 
crack depth and the height of the element (s=alh), defined as 
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The flexibility coefficient of the intact element can be derived 

as 
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The additional flexibility coefficients are 
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The flexibility matrix of the intact element is 

Similarly the coefficients C„ ' can be expressed in matrix 
form as 
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The total flexibility coefficients for the element with crack are 
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The total flexibility matrix for the element with crack can be 
expressed as 
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From the equilibrium conditions the following relationship 
holds 
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Where 

"-I -L 1 0" 

0-101 
The stiffness matrix of the intact element can be written as 
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The stiffness matrix of the cracked element can be written as 



(15) 



The Euler-Bernoulli's beam element is used in the formulation. 
It is assumed that the crack present in the element affect only 
stiffness but not the mass. The cracked beam element is 
introduced with concept of linear fracture mechanics. The free 
vibration analysis of the beam is performed using MATLAB 
code. The consistent mass matrix and stiffness matrix for beam 
element is used, to analyze the intact beam. The algorithm 
proposed in the paper follows. 

III. Algorithm For Damge Detection 
The proposed algorithm for damage detection is 

> Forward problem 

1 . Model the structure for solving eigenvalue problem in finite 
element method. 

2. Use model updating to match the responses from finite 
element method and actual response. 

3. Obtain Normalized mode shapes and frequencies from the 
eigen solution. Curvature of mode shapes is obtained and 
used for damage location. 

Damage response surface function, Error damage response 
surface function and Resultant damage response function 
are established as a function of location of crack from fixed 
end (Lcr) and crack depth ratio (a/h). Few typical fit can be 
mathematically written as 

> Quadratic fit surface polynomial 
Curd = A+ BL rp + CY+ DL 2 +EY 2 +FL rB Y 

CK CK LK 

& 
Error=A 1+ B 1 L CR+ C 1 Y + D 1 L CR 2 + E i r- + F 1 L CR Y 

> Cubic fit surface polynomial 

Curd = A+BL CR +C Y+DL CR 2 +EY 2 +FL CR Y 

+GL tR 3 +m- + IL CR 2 Y + 'L CR Y 2 

& 

Error Curd = A, + B, L rR + C , Y + D,L 2 E, Y 2 

1 1 Uv 1 1 L.K 1 

+ F 1 L CR Y + G 1 L CR VH 1 V + I 1 L CR 2 Y + ' 1 L CR Y 2 

Where Curd = curvature. Y= (a/h) = Crack depth ratio & L = 

Crack location from fixed end 

The resultant curvature and error curvature are calculated as 

follows 
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Resultant Curd = Curd - error. 
Error Curd = Curd new - Curd 

> Inverse problem 

To simulate real life situation Noise added in mode shapes. 
Crack location is located by wavelet transform or curvature 
of mode shape. Crack depth ratio (a/h) is expressed as a 
function of location of crack from fixed end (Lcr) and Curvature 

i.e. 
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Construction and plotting of Damage response surface 
function, Error damage response surface function and 
Resultant damage response function. .Percentage error is 
established from the numerical simulated data. 

IV. NUMERICALSIMULATION 

To implement the proposed algorithm, a cantilever beam is 
considered with following geometrical parameters. The beam is 
0.3 m long with cross section 0.0225 m x0.013 m respectively.. 
The material properties taken are: modulus of elasticity 175 
GPa, mass density 7.8 kg/m3, and Poisson ratio 0.3. The crack 
depth 'a' and crack location L are varied during the analysis. 
Based on the proposed algorithm a code was generated. 
Numerical simulation obtained by finite element by 
introducing crack at a location with variation in crack depth. 
The crack locations are varied from to 300 in a step of 50 
while the crack depth ratio varied from intact to 0.6. The code 
prompts user to solve the problem with or without Guyan 
reduction technique. The eigenvalue problem is solved to 
obtain modal parameters. 
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Figure 1. Curvature with varying Crack depth (for length 100 and 
150 (mode 1)) 

The curvature mode shapes are evaluated. The variation 
of curvature increases as the crack depth increases. 
Additional data is generated by meshing the numerical data 
with spline function. 
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Figure 2. Spline fitted Curvature Surface plot. 
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Figure 3. Damage curvature response surface plot 

The damage surface for the curvature is generated using 
response surface method. MATLAB code is used for 
obtaining coefficients of quadratic and cubic polynomial fit. 
Such polynomial fit is obtained for first five mode shape. 
Equation below is a quadratic polynomial fit for mode 3. 

Curd=1.23e-03+5.51e-05*L CR +7.78e-03*Y 

+2.02e-07*L CR 2 +4.27e-2*Y 2 +1.47e-05*L CR Y 

Construction and plotting of Damage response surface 
function, Error damage response surface function and 
Resultant damage response function as a function of location 
of crack from fixed end (Lcr) and crack depth ratio (a/h) from 
numerical data. Similarly, Construction and plotting of 
Damage response surface function, Error damage response 
surface function and Resultant damage response function. 
Crack depth ratio (a/h) is expressed as a function of location 
of crack from fixed end (Lcr) and Curvature i.e. 

3. 

— = F[ (Lcr, Curvature] where : =1.2.3.... X 

The inverse problem is used for damage detection. To 
simulate real life situation Noise added in Normalize mode 
shapes. The wavelet transform (figure 4) is used for locating 
crack the location of crack is verified by plotting the curvature 
mode shape. Even multiple cracks are also located by Wavelet 
transform. 
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Figure 4. Error damage curvature response Surface plot 




Figure 5. Crack detection at element number 75 (with crack depth 
0.05) 

Resultant damage response polynomial function in terms of 
the Crack depth ratio (a/h) is expressed as a function of 
location of crack from fixed end (Lcr) and Curvature from 
numerical data. This function is used to establish Crack depth 
ratio. 

TABLE 1 
PERCENTAGE ERROR IN PREDICTING CRACK DEPTH RATIO 
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Table 1 shows Percentage error in predicting Crack depth 
ratio for few cases. 



The percentage error indicates the lower percentage error for 
most of the crack depth ratio except for the lower crack depth 
ratio. The percentage error is reduced in case of cubic fit but 
error increases at lower crack depth. 

iv. Conclusions 

The present study addresses on the identification of the 
presence, location and size of a crack in a cantilever beam 
structure by extracting mode shapes and frequencies from 
numerical simulation. The eigenvalue problem is solved to 
obtain the modal parameters. The algorithm presented based 
on curvature mode shape which is used for damage detection. 
The new concept of error surface plot is introduced for 
efficient and accurate prediction of crack location and size. 
The wavelet transform shows the sharp peaks at the crack 
location and depth is established by response surface fitting 
function. The error in the prediction of crack depth ratio is 
higher for lower crack depth ratio. The technique is applicable 
to any linear structure that can be accurately modeled using 
finite element method. 
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